Analytic results for the three-sphere swimmer at low Reynolds number 
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The simple model of a low Reynolds number swimmer made from three spheres that are connected 
by two arms is considered in its general form and analyzed. The swimming velocity, force-velocity 
response, power consumption, and efficiency of the swimmer are calculated both for general defor- 
mations and also for specific model prescriptions. The role of noise and coherence in the stroke cycle 
is also discussed. 
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I. INTRODUCTION 

There is a significant compHcation in designing swim- 
mers at small scale as they have to undergo non- 
reciprocal deformations to break the time-reversal sym- 
metry and achieve propulsion at low Reynolds number 
[H. While it is not so difficult to imagine constructing 
motion cycles with the desired property when we have a 
large number of degrees of freedom at hand — like nature 
does, for example — this will prove nontrivial when we 
want to design something with only a few degrees of free- 
dom and strike a balance between simplicity and func- 
tionality, like most human-engineered devices Re- 
cently, there has been an increased interest in such de- 
signs 
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an interesting example of such robotic micro-swimmers 
has been realized experimentally using magnetic colloids 
attached by DNA-linkers While constructing small 
swimmers that generate surface distortions is a natural 
choice, it is also possible to take advantage of the gen- 
eral class of phoretic phenomena to achieve locomotion, 
as they become predominant at small scales [Toj . 

Here we consider a recently introduced model for a 
simple low Reynolds number swimmer that is made of 
three linked spheres and present a detailed analy- 
sis of it motion. Unlike the Purcell swimmer [2] that 
is difficult to analyze because it takes advantage of the 
rotational degrees of freedom of finite rods that move 
near each other , the three-sphere swimmer model is 
amenable to analytical analysis as it involves the trans- 
lational degrees of freedom in one dimension only, which 
simplifies the tensorial structure of the fluid motion. We 
present closed form expressions for the swimming veloc- 
ity with arbitrary swimming deformation cycles, and also 
use a perturbation scheme to simplify the results so that 
the study can be taken further. We examine various 
mechanical aspects of the motion including the pattern 
of the internal forces during the swimming, the force- 
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velocity response of the swimmer due to external loads, 
the power consumption rate, and the hydrodynamic effi- 
ciency of the swimmer. Finally, we consider the role of 
the phase difference between the motion of the two parts 
of the swimmer and propose a mechanism to build in a 
constant (coherent) phase difference in a system that is 
triggered from the outside. We also discuss the effect of 
noise on the swimming velocity of the model system. We 
also note that the three-sphere low Reynolds swimmer 
has been recently generalized to the case of a swimmer 
with a macroscopic cargo container [2l| . and a swimmer 
whose deformations are driven by stochastic random con- 
figurational transitions [l^l . 

The rest of the paper is organized as follows. In Sec. 
|TT]the three-sphere low Reynolds number swimmer model 
is introduced in a general form and a simplified analysis 
of its swimming is presented. This is followed by a de- 
tailed discussion of its swimming velocity in Sec. IIIIl 
with a focus on a few particular examples of swimming 
stroke cycles. Section |IV] is devoted to a discussion on 
internal stresses and forces acting during the swimming 
cycle, and Sec. |V] studies the force-velocity response of 
the swimmer. The power consumption and efficiency of 
the model swimmer are discussed in Sec. I VII followed by 
concluding remarks in Sec. IVIII Appendix [X] contains 
the closed form expression for the swimming velocity of 
the general asymmetric swimmer, which is the basis of 
some of the results discussed in the paper. 



II. THREE-SPHERE SWIMMER: SIMPLIFIED 
ANALYSIS 



We begin with a simplified model geometry which con- 
sists of three spheres of radii (i = 1, 2, 3) that are sep- 
arated by two arms of lengths Li and L2 as depicted in 
Fig. [TJ Each sphere exerts a force fi on (and experi- 
ences a force —fi from) the fiuid that we assume to be 
along the swimmer axis. In the limit Ui/Lj <C 1, we can 
use the Oseen tensor [2^, to relate the forces and the 
velocities as 
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Note that in this simple one dimensional case, the ten- 
sorial structure of the hydrodynamic Green's function 
(Oseen tensor) does not enter the calculations as all the 
forces and velocities are parallel to each other and to the 
position vectors. The swimming velocity of the whole 
object is the mean translational velocity, namely 



(4) 



We are seeking to study autonomous net swimming, 
which requires the whole system to be force-free (i.e. 
there are no external forces acting on the spheres). This 
means that the above equations are subject to the con- 
straint 



h+f2+h= 0. 



(5) 



Eliminating /2 using Eq. ([5]), we can calculate the swim- 
ming velocity from Eqs. ([T]), (g]), ([3]), and (g]) as 
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where the subscript denotes the force-free condition. To 
close the system of equations, we should either prescribe 
the forces (stresses) acting across each linker, or alterna- 
tively the opening and closing motion of each arm as a 
function of time. We choose to prescribe the motion of 
the arms connecting the three spheres, and assume that 
the velocities 



Li = V2 

L2 = V3 
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are known functions. We then use Eqs. ((T]), and 
(ini) to solve for fi and /a as a function of Li and L2- 
Putting the resulting expressions for fi and /s back in 
Eq. ([5]), and keeping only terms in the leading order in 
ai/Lj consistent with our original scheme, we find the 
average swimming velocity to the leading order. 




FIG. 1: Schematics of the three-sphere swimmer. The two 
arms can open and close in a prescribed form, and this could 
lead to locomotion if the swimming strokes are nonreciprocal. 



III. SWIMMING VELOCITY 

The result of the above calculations is the lengthy ex- 
pression of Eq. (jAip reported in Appendix El This re- 
sult is suitable for numerical studies of swimming cy- 
cles with arbitrarily large deformations. For the simple 
case where all the spheres have the same radii, namely 
a = ai = a2 = 0,3, Eq. ^ simplifies to 




(9) 



plus terms that average to zero over a full swimming 
cycle. Equation [5] is also valid for arbitrarily large defor- 
mations. 

We can also consider relatively small deformations and 
perform an expansion of the swimming velocity to the 
leading order. Using 



Li = £1+ ui, 

L2 = £2+ U2, 



(10) 
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in Eq. (|A1[) . and expanding to the leading order in Ui/£j, 
we find the average swimming velocity as 
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In the above result, the averaging is performed by time 
integration in a full cycle. Note that terms proportional 
to u\U\, U2U2, and ^1-02 -t- U1U2 are eliminated because 
they are full time derivatives and they average out to zero 
in a cycle. Equation ([T2|) clearly shows that the average 
swimming velocity is proportional to the enclosed area 
that is swept in a full cycle in the configuration space [i.e. 
in the (ui,U2) space]. This result, which is valid within 
the perturbation theory, is inherently related to the ge- 
ometrical structure of theory the low Reynolds number 
swimming studied by Shapere and Wilczek [2^. Natu- 
rally, the swimmer can achieve higher velocities if it can 
maximize this area by introducing sufficient phase dif- 
ference between the two deformation cycles (see below). 
We also note that the above result is more general than 
what was previously considered in Ref. Q, which cor- 
responded to the class of configurational changes that 
happen one at a time, i.e. spanning rectangular areas in 
the configuration space [25j . 

We can actually obtain Eq. from a rather gen- 

eral argument. Since the deformation of the arms is pre- 
scribed, the instantaneous net displacement velocity of 
the swimmer should take on a series expansion form of 
v{t) = Aiiii + BijiiiUj + CijkUiUjUk + • ■ ■ , where the co- 
efficients Ai, Bij, Cijk, etc. are purely geometrical pref- 
actors {i.e. involving only the length scales o^'s and ^i's). 



3 



Terms of higher order than one in velocity will have to be 
excluded on the grounds that in Stokes hydrodynamics 
forces are linearly dependent on (prescribed) velocities, 
and then the velocity any where else is also linearly pro- 
portional to the forces, which renders an overall linear 
dependency of the swimming velocity on set velocities. 
Moreover, higher order terms in velocity would require a 
time scale such as the period of the motion to balance 
the dimensions, which is not a quantity that is known 
to the system at any instant (i.e. would require nonlo- 
cal effects). Since the motion is periodic and we should 
average over one complete cycle to find the net swim- 
ming velocity, we can note that the only combination 
that survives the averaging process up to the second or- 
der is U1U2 — U1U2, which yields Eq. (fT2|) . Note that this 
argument works even if the spheres have finite large radii 
that are not small comparable to the average length of 
the arms. 

It is instructive at this point to examine a few explicit 
examples of swimming cycles for the three-sphere swim- 
mer. 



A. Harmonic Deformations 

Let us consider harmonic deformations for the two 
arms, with identical frequencies uj and a mismatch in 
phases, namely 



u 



',i{t) — di cos{ujt + ipi), 
U2{t) — d2 cos((jji + ^2)- 



(14) 
(15) 



The average swimming velocity from Eq. (|12p reads 



K 



did2L0 8\-a{ipi - 1P2). 



(16) 



This result shows that the maximum velocity is obtained 
when the phase difference is 7r/2, which supports the pic- 
ture of maximizing the area covered by the trajectory of 
the swimming cycle in the parameter space of the defor- 
mations. A phase difference of or tt, for example, will 
create closed trajectories with zero area, or just lines. 



B. Simultaneous Switching and Asymmetric 
Relaxation 

Thinking about practical aspects of implementing such 
swimming strokes in real systems, it might appear dif- 
ficult to incorporate a phase difference in the motion of 
the two parts of a swimmer. In particular, for small scale 
swimmers we would not have direct mechanical access to 
the different parts of the system and the deformations 
would be more easily triggered externally by some kind 
of generic interaction with the system, such as shining 
laser pulses. In this case, we need to incorporate a net 
phase difference in the response of the two parts of the 
system to simultaneous triggers. This can be achieved 



if the two parts of the system have different relaxation 
times. To illustrate this, imagine that the arms of the 
swimmer could switch their lengths via exponential re- 
laxation between two values oi £i — di/2 and £i + di/2 
back and forth as a switch is turned on and off. The 
deformation can be written as 
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for i = 1,2, where r^'s are the corresponding relaxation 
times, and T is the (common) period of the switchings. 
We find 
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The above function is a smooth and monotonically decay- 
ing function of both x and y that is always positive, and it 
has the asymptotic limits ^(0, 0) = |, ^{x, y — + 00) — 0, 
and J-{x — !■ 00, y) = 0. Here, the phase mismatch is 
materialized in the difference in the relaxation times, de- 
spite the fact that the deformations are switched on and 
off simultaneously. 



C. Noisy Deformations 

Another important issue in practical situations is the 
inevitability of random or stochastic behavior of the de- 
formations. In small scales, Brownian agitations of the 
environment become a significant issue, and we would like 
to know how feasible it is to extract a net coordinated 
swimming motion from a set of two noisy deformation 
patterns. Using the Fourier transform of the deforma- 
tions {i — 1, 2) 
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we can calculate the time-averaged swimming velocity of 
Eq. (ini) as 

^ K I f doj 
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(21) 

For deformations that have discrete spectra [i = 1,2), 
namely 



i{t) = ^ din COs{uJnt + ipin), 



(22) 



we find 



V^o = y ^ di„d2„w„ sin((^i„ -(p2n)- (23) 
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FIG. 2: (color online.) The dimensionless forces exerted on 
the fluid at the locations of the spheres [defined in Eq. (|30|l ] 
as functions of the relative size of the two arms. Dashed line 
(red) corresponds to Sphere 1, dotted line (blue) corresponds 
to Sphere 2, and solid line (green) corresponds to Sphere 3. 



This shows that a net swimming is the result of coor- 
dinated motions of the different modes in the frequency 
spectrum, and the net velocity is the sum of the individ- 
ual contributions of the different modes. As long as we 
can achieve a certain degree of coherence in a number of 
selected frequencies, we can have a net swimming despite 
the noisy nature of the deformations. 



IV. INTERNAL FORCES AND STRESSES 



The interaction of the spheres and the medium involves 
forces as these parts of the swimmer make their ways 
through the viscous fluid when performing the swimming 
strokes. We have calculated these forces for a general 
swimmer, but since their expressions are quite lengthy, 
we choose to present them only in the particular case 
where all the spheres have equal radii, namely ai = a2 — 
03 = a. For arbitrarily large deformations, we find 
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up to terms that average to zero over a full cycle. One 
can check that the above expressions manifestly add up 
to zero, as they should. For small relative deformations, 
we find the following expressions for the average forces 



(a) 



(b) 



Swimming velocity 



FIG. 3: (color online.) The distribution of the average forces 
(denoted by arrows) exerted on the different spheres depend- 
ing on the relative size of the arms. Case (a) corresponds to 
i\ = I2, case (b) corresponds to ii <C ^2, and case (c) corre- 
sponds to l\ S> ^2. The middle sphere always experiences a 
force that pushes it in the direction of swimming irrespective 
of the structure of the swimmer. Note that the above different 
cases do not correspond to the instantaneous forces exerted 
on each sphere during the swimming cycle. 
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These forces are all proportional to the net average swim- 
ming velocity, with proportionality constants that de- 
pend on li and £2- Therefore, we can write a generic 
form 



(30) 



for the forces in terms of the dimensionless factors gi. 
These dimensionless forces are plotted in Fig. [3] as func- 
tions of the ratio between the lengths of the two arms, 
which is a measure of the asymmetry in the structure of 
the swimmer. 

It is instructive to examine the limiting behaviors of 
the forces as a function of the asymmetry. For £1 <g; £2 
we have gi — — 1 , g2 — — 1 , and g^ — 2 which means that 
the two closer spheres are in line with each other and 
the third sphere that is farther apart exerts the opposite 
force. The same trend is seen in the opposite limit £1 ^ 
H.2, where we have = 2, (72 = — 1, and g^ = —1. In the 
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symmetric case where ii = £2 we have ffi = ffa = f and 
(72 = Note that 32 is always negative, which means 

that the middle sphere always experiences a force from 
the fluid that pushes it in the direction of swimming. 
The distribution of forces exerted on the spheres in these 
three different limits is shown schematically in Fig. [31 

It is interesting to note that in the symmetric case 
the forces are distributed so that their net dipole mo- 
ment vanishes and the first non- vanishing moment of the 
forces becomes the quadrapolc moment. This will cause 
the net fall-off of the velocity profile at large separations 
r to change from l/r^ (force dipoles) into (force 
quadrapoles) . This behavior can be explained by a more 
general symmetry argument p^ . 



V. FORCE-VELOCITY RELATION AND STALL 
FORCE 

The effect of an external force or load on the efficiency 
of the swimmer can be easily studied within the linear 
theory of Stokes hydrodynamics. When the swimmer is 
under the effect of an applied external force F, Eq. ^ 
should be changed as 



/l + /2 + /3 = F. 



(31) 



Following through the calculations of Sec. |TT] above, we 
find that the following changes take place in Eqs. ((T|), 
(HD, ©, and Q: 
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in Eq. (jAip . which together with correction coming from 
Eq. ([35)) leads to the average swimming velocity 



V{F) = Vo + 



F 



(39) 



to the leading order, where aji is an effective (renormal- 
ized) hydrodynamic radius for the three-sphere swimmer. 
To the zeroth order, we have = ^(oi + a2 + 03) for 



the general case and there are a large number of correc- 
tion terms at higher orders that we should keep in order 
to be consistent in our perturbation theory. Instead of 
reporting the lengthy expression for the general case, we 
provide the expression for oi =02 = 03 = a, which reads 
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We can also expand the deformation up to second order 
in Ui and average the resulting expression. The result 
of this calculation is also lengthy and not particularly 
instructive, and is hence not reported here. 

The force- velocity relation given in Eq. (139]) . which 
could have been expected based on linearity of hydrody- 
namics, yields a stall force 
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(41) 



Using the zeroth order expression for the hydrodynamic 
radius, one can see that this is equal to the Stokes force 
exerted on the three spheres moving with a velocity Vq. 



VI. POWER CONSUMPTION AND 
EFFICIENCY 

Because we know the instantaneous values for the ve- 
locities and the forces, we can easily calculate the power 
consumption in the motion of the spheres through the 
viscous fiuid. The rate of power consumption at any time 
is given as 

V = hvi + f2V2 + hv3 = /i(-ii) + 73(^2), (42) 

where the second expression is the result of enforcing the 
force- free constrain of Eq. Using the expressions for 
/i and /3 as a function of Li and L2, we find 
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(43) 



for ai = a2 — a^ — a. 

We can now define a Lighthill hydrodynamic efficiency 

as 
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for which we find to the leading order 
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where d - ^ + f,-hf; + T;fr2^ C'a - 1-^17 + ^ + 17717' 
and C3 = 1 — ^^ — 5^ + 1 17^17 ■ It is interesting to note 
that for harmonic deformations (with single frequency) 
Eq. (|l5)) is independent of the frequency and scales like 
a?<P /l'^, which reflects the generally low hydrodynamic 
efficiency of low Reynolds number swimmers. In this 
case, it is possible to find an optimal value for the phase 
difference that maximizes the efficiency @. 

VII. CONCLUDING REMARKS 

We have considered the simple model of a low Reynolds 
number swimmer, which is composed of three spheres 
that are linked by two phantom arms with negligible hy- 
drodynamic interaction. Assuming arbitrary prescribed 
motion of the two arms, we have analyzed the motion 
of the swimmer and provided explicit expressions for the 
swimming velocity and other physical characteristics of 
the motion. 

The simplicity of the model allows us to study the 
properties of the swimmer in considerable details using 
analytical calculations. This is a great advantage, as it 
can allow us to easily consider complicated problems in- 
volving such swimmers and could hopefully lead to new 
insights in the field of low Reynolds number locomotion. 



An example of such studies has already been performed 
by Pooley et al. who considered the hydrodynamic in- 
teraction of two such swimmers and found a rich vari- 
ety of behaviors as a function of relative positioning of 
the two swimmers and their phase coherence iMj. This 
can be further generalized into a multi-swimmer system, 
and the collective floe behavior of such systems can then 
be studied using a "realistic" model for self-propellers at 
low Reynolds number that is faithful to the rules of the 
game. Knowing something about the internal structure 
of a swimmer will also allow us to study the synchroniza- 
tion problem more systematically 
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APPENDIX A: SWIMMING VELOCITY FOR 
ARBITRARY DEFORMATIONS 

For the general case of a three-sphere swimmer based 
on the schematics in Fig. [Jl we obtain the average swim- 
ming velocity to the leading order as 
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This expression can be used in numerical studies of the geometric characteristics, 
swimming motion for arbitrarily large deformations and 
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